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Abstract 

The Bose-Hubbard Hamiltonian of spin-2 cold bosons with repulsive interaction in an optical lattice 
is proposed. After neglecting the hopping term, the site-independent Hamiltonian and its energy 
eigenvalues and eigenstates are obtained. We consider the hopping term as a perturbation to do the 
calculations in second order and draw the phase diagrams for different cases. The phase diagrams 
show that there is a phase transition from Mott insulator with integer number bosons to superfluid 
when the ratio co/t (cq is the spin- independent on-site interaction and t the hopping matrix element 
between adjacent lattice sites) is decreased to a critical value and that there is different phase boundary 
between superfluid and Mott insulator phase for different Zeeman level component in some ground 
states. We find that the position of phase boundary for different Zeeman level component is related 
to its average population in the Mott ground state. 
PACS number(s):03.75.Kk, 03.75.Lm, 03.75.Mn, 32.80.Pj 



I. INTRODUCTION 



During recent years, the superfluid-Mott insulator transition becomes a significant topic of 
atomicphysics, which has been extensively studied in the context of 4 He absorbed in the porous 
media^]. The detailed study of superfluid-Mott insulator transition becomes possible when 
quantum phase transition from a superfluid to a Mott insulating ground state was observed in 
a Bose-Einstein condensate stored in a three dimensional optical lattice potential by Greiner 
et al.0. The optical lattices used to confine cold atoms are induced by the ac Stark effect of 
interfering laser beams. The dynamics of the bosonic atoms in optical lattices realizes a Bose- 
Hubbard model which predicts phase transition from a superfluid(SF) to Mott insulator(MI) 
at low temperature with increasing the ratio of the on site interaction Co to the hopping matrix 
element t. 

Besides many experimental efforts made to realize superfluid-Mott insulator transition, a 
la rg e nu mta of theo.et.ca> >i_ appeared fl Q fl ^ Q. Ref.fl mad e an appropriate 
mean-field approximation to the Hamiltonian of spinless or polarized bosons in an optical lattice 
to describe the zero-temperature phase transition from the superfluid to the Mott-insulating 
phase analytically and calculated the phase diagram. In this paper a first-order approximation 
to the dispersion of the density fluctuations shows that the system indeed goes from a gapped 
to a gapless phase. 

Since trapping a Bose condensate by purely optical means is realized several years ago[ 
which liberates the internal degrees of freedom of spinor Bosons that are frozen in a magnetic 
trap, much attention is attracted on the study of spinor Bose condensate, which has multiple 
components. Hoj^ and Ohmi and MachidajlQ] have researched a spin-1 BEC and obtained the 
general theoretical frameworks respectively. Ciobanu,Yip and Ho[ljJ have studied an F = 2 
spinor bose condenate and drawn the phase diagrams. Koashi and Ueda have obtained the exact 
eigenstates of spin-1 and spin-2 vertical BEC and discussed their magnetic response 3, Q|- 

The next natural step after the detailed research of spinor bose condensates with free internal 
freedom is the study of the spinor multi-component BEC in a optical lattice — a new background. 
Demler and Zhou have studied spin-1 Bose atoms in an optical lattice and obtained several 
unique properties Ref.j^J studied the spin-1 bosons interacting antiferromagnetically in an 
optical lattice with the mean-field approximation method and obtained the phase diagram 
showing a transition from Mott insulator to superfluid, which is signed by the appearing of the 
superfluid order parameter i.e. the density fluctuations. But the spin-2 cold Bose atoms in an 
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optical lattice have not been researched as yet, which is the principal task of this paper. 

In this paper, we study the transition from Mott insulator to superffuid of spin-2 Bose atoms 
with repulsive interaction in terms of mean-field approximation method. First, ignoring the 
hopping term of the Hamiltonian, we get the eigenstates and eigenvalues and find the ground 
states of the system for different cases. Then, we consider the hopping term as a perturbation 
to do calculations in second order and draw the phase diagrams. 

II. THE MODEL 

We consider bosons with hyperfine spin F = 2, such as 23 Na, 87 Rb or 85 Rb, in an optical 
lattice. Adding the external periodic potential term V(r) = Vo(sin 2 fca; + sin 2 ky + sin 2 kz), 
where Vq is a tunable amplitude and k the wave vector of the laser light, to the Hamiltonian 
of spin-2 bosons we get the Hamiltonian of spin-2 bosons in an optical lattice in second 
quantized form as follows 

+ | £(*£(*k)<tf*/j) 2 + c 2 ^l(2a; 2a'|00)(00|2/3; 2/7)¥ /3 *>] (1) 

where m is the atomic mass, \l/+2, ^-2 are the five-component field operators corresponding 
to the sublevels m F = +2, — 2 of the hyperfine state F = 2, \i is the chemical potential, c , 
c\ and 62 are related to s-wave scattering length 0-0,02 and 04 of the two colliding bosons with 
total angular momenta 0, 2 and 4 by cq = 47rfr 2 (3a4 + 4a2)/7m, c\ = 4irh 2 (a^ — 0,2) /7m and 62 = 
47rft 2 (3o 4 - 10o 2 + 7o ) /7m, < 2a; 2a'|00 > and < 00|2/3; 2(5' > are Clebsch-Gordan coefficients. 
F a (a = x, y, z) are 5x5 spin matrices satisfying the commutation relation [F a , Fp] = Eap^F^. 

For a single atom in the periodic potential, the energy eigenstates are Bloch states. In the 
tight-binding limit, we can superpose the Bloch states to get a set of Wannier functions, which 
are localized on an individual lattice site. Expanding the field operators in the Wannier basis 
and keeping only the lowest vibrational states, ^f a = J2ibi a w(r — r^), Eq.((H) reduces to the 
Bose-Hubbard Hamiltonian 

H = -tJ2 &I&i-A*£A* + ?£*(*-i) + T E^-^ + yEM- (2) 

<i,j> i Z i Z i i 

where F; = bl a F af3 b if3} and hi = Y. a b\ a bi a . t = - J drw*(r)(-h 2 V 2 /2m + V(r))wj(r) is the 
hopping matrix element between adjacent sites % and j. Ci is on-site inter-atom interaction 
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defined by q = q / dr|i0j(r)| 4 , where the Hubbard approximation has been used to approximate 
the multi-center integral as a single-center one. Si + = Sj_ = (&j ) 2 /2 — &li&J_i + b\ 2 bl_ 2 , when 
applied to the vacuum, creates a spin-singlet "pairs". The properties of the operators will 
be discussed in detail in the following section. The fourth and fifth terms in the equation are 
adding terms comparing with the Bose-Hubbard Hamiltonian of spinless or polarized bosons [3] . 
By contract with that of spin-1 bosons jjj], the fifth term is an adding term and 6fii in the fourth 
term replaces 2nj. 

In the limit co/t — ► oo, the hopping term in the Hamiltonian can be neglected, so the 
Hamiltonian is reduced to a diagonal matrix with respect to sites. Omitting the site index, the 
singe-site Hamiltonian is 

ho = -fin + ^-h(n - 1) + ^-(F 2 - Qh) + (3) 

Zi Zi 

whose eigenstates and eigenvalues will be discussed in following section. 

To study the quantum transition, we use the mean-field approximation to the hopping 
term and consider it as a perturbation. Here we introduce the superfluid order parameter 
(p a =< b ia >= y/n s f( a , where n s f is the superfluid density and ( a is a normalized spinor 
CaCa = 1- The hopping term is decoupled as b\ a bj a « {4> a b\ a + <p a bj a ) — (fi a 4> a - So the total 
hopping term becomes the product of a site-independent term and the total number of the 
sites. We can consider only a single site because the Hamiltonian of every site is identical 
in the homogenous case. In the single site, the resulting mean-field version of the hopping 
Hamiltonian can be written as 

h = zt((j) Q b a + <f>* a b a - <f) a (f) a ) (4) 

where z is the number of the nearest-neighbor sites and the site-index is omitted. So the 
Hamiltonian of a single site is 

h = h + h 1 (5) 
When the ratio c /t is very large, we can consider h± as a perturbative term. 

III. THE ENERGY EIGENVALUES AND THE GROUND STATES 

Before perturbative calculations we should get the eigenvalues and eigenstates of Eq.Q, i.e., 
to sovle the equation 

h ^jj = e <°ty (6) 
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In Eq.flSJ), the operaters SU and 5"_ satisfy the SU(1, 1) commutation relations together with 
S z = (2n + 5)/4, namelyflQ, 

[S Z ,S ± ]=±S ± , {S + ,S-] = -2S z (7) 

~ 2 J 

and the Casimir operator S commuting with and reads 

S 2 = + 5 2 2 - 4 (8) 
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The eigenvalues of the mutual eigenstates for S and S z are {S(S — 1), S z } with 5 = (2n + 
5)/4 (n = 0,1,2,...) and S z = S + n s (n s = 0,1,2,...), which guarantees that S + S- = 
Si — S z — S 2 is positive semidefinite. The new quantum numbers n s and no are introduced and 
the operator S + raises n s by one while the relation n = 2n s + uq holds, where n is the total 
number of bosons in a single site. 

The operators S± commute with the spin operator F and the magnetic quantum number 
operator F z , So the energy eigenstates can be classified according to quantum numbers no and 
n s , total spin F, and magnetic quantum number F z . Hence, the eigenstates if) are denoted as 
\no, n s , F, F z \ A > where A labels orthonormal degenerate states. The energy eigenvalue is given 
by 

e (o) = _ un + ^ n ( n _ i) + C _L\F(F + i) _ 6n ] + ^n s (n - n, + |) (9) 

where the relation n = 2n s +n is used. So Mott states can be expressed as Hi \n>o, n s , F, F z \ A >j, 
where % is lattice site index. For homogenous case, the zeroth order total energy is 
J2i £^ = Niey\ where Ni is the number of lattice sites. 

From Ref. |l3l| we know that the energy eigenstates \uq, n s , F, F z \ A > can be represented as 

AF( J(2)t\n 2 o p. , (U\\ni2( A^\ n 22( /j( 3 )t\n 30 /' A(3)t 



(F_)^(^ T )"-P (n3=0) (6D ni2 (^2 T 22 (Ar) n30 (Ar T"\vac> (10) 
where 

4 2)t = ^[( & S) 2 " 2^ f -i + 26 2 6l 2 ] (11) 

4 2)t = -L[2V2bibi - V3(b{) 2 ] (12) 



i 3)t = -±={V2(bl) 3 - 3V2b\blbU + 3V3(b\)%U + 3V36 2 (6_0 2 - 6v^&S&l 2 ] (13) 
4 3)t = ^I) 3 - VQbW + 2(6 2 )V_ 1 ] (14) 
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P(n s =o) is the projection onto the subspace with n s = 0; n 12 , n 20 , ^22, ^30 = 0, 1, 2, 00, n 33 = 
0, 1, and AF = 0, 1, 2F that are related to {n , n s , F, F z } as 



n = "12 + 2n 22 + 3n 30 + 3n 33 



(15) 



n s = n 20 



(16) 



F 



2n i2 + 2n 22 + 3n : 



33 



(17) 



F z 



F-AF 



(18) 



To find the ground states of Eq.fjHJ), we minimize its energy eigenvalue, i.e., Eq.®. We note 
that there isn't magnetic quantum number F z in Eq.fjHJ), which means the eigenstates with 
different magnetic quantum number but the same other quantum numbers are degenerate. For 
simplicity, in this paper we only investigate the states with the highest magnetic quantum 
number, i.e. the states with F z = F. when n — 1, it is obvious that the ground state is 
|1, 0, 2, 2; A >, which means only one single atom with the highest magnetic quantum number 
for the state. For n > 2, it is more complicated, which will be discussed in detail as follow: 

1. Ci > 0, c 2 > 0. In this case, if n s = and F = 0, the energy eigenvalue has its minimum. 
But that is not always the case because there are some forbidden values [12I Q|, i.e., 
F — 1,2, 5, 2tiq — 1 are not allowed when no = 3k(k £ Z), and F — 0,1, 3, 2uq — 1 are 
forbidden when no = 3k ± l(k £ Z), according to which, it is classified into three cases. 

(a) For n = 3k(k £ Z), \n, 0, 0, 0; A > with n s = and F = simultaneously, which is 
not within forbidden values, is the ground state 

(b) For n = 3k — l(k £ Z), the state with n s = and F = simultaneously is not 
allowed. If n s = 0, the lowest allowed value of total spin F is 2. F = is not 
forbidden when no is 3k(k £ Z). n s = 1, due to the relation n = 2n s + no, is the 
lowest value satisfying the condition. So the state with n s = 1 and F = is a possible 
ground state. There is competition between the third term and the fourth term in 
Eq.(jnj),i.e., competition between contributions of total spin and singlet "pairs" to 
eigenenergy. Comparing both eigenenergies for the two cases F = 0, n s = 1 and 
F = 2,n s = 0, we get the ground state, (i)|n — 2, 1,0,0; A > for c\ > c 2 (2n + 1)/15 
and (ii) \n, 0, 2, 2; A > for c x < c 2 (2n + 1) /15. 

(c) For n = 3k + l(k £ Z), it is similar to last case that there is competition between 
contributions of total spin and the singlet "pair" to eigenenergy in this case. But 
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there is a difference between the two cases in that if F = 0, n s is at least 2. Therefore, 
the ground state is (i)|n - 4, 0, 2, 2; A > for ci > 2c 2 {2n - 1)/15; (ii) |n, 0, 2, 2; A > 
for C\ < 2c 2 (2n — 1)/15, which is obtained by comparing the eigenenergy of these 
states. 

2. ci > 0, c 2 < 0. In this case, the eigenenergy is the lowest when F is at the lowest value 
and n s at the highest value. There is a state with F = and n s = n/2 when n is even, 
but there is not when n is odd. So, when n is different, there are two cases as follow: 

(a) When n is even, the ground state is \n, n/2, 0, 0; A > 

(b) When n is odd, n s has the highest value (n — l)/2. But F is not zero when n s = 
(n — l)/2. Alternatively, there is another case that F = and n s = (n — 3)/2. In 
the two cases, the state that has lower eigenenergy is the ground state. Hence, the 
ground state is (i) |1, (n - l)/2, 2, 2; A > for c x < 7|c 2 |/15; (ii) |3, (n — 3)/2, 0, 0; A > 
for ci > 7|c 2 |/15. 

3. For Ci < 0, c 2 > 0, the eigenstate with F = 2n and n s = 0, which has the lowest 
eigenenergy, is the ground state. So the ground state is \n, 0, 2n, 2n; A >. 

4. For c\ < 0, c 2 < 0, to get the ground state, we minimize Eq.Q. If we skip over the 
fact for the moment that the singlet "pair" number n s is an integer, then the condition 
to minimize the energy function is given by 

_ 10ci(4rc+ 1) -c 2 (2n + 3) 
Us 80ci-4c 2 [ ' 

because of the fact that the singlet "pair" number must be an integer, we write n s in 
terms of the closet integer number n° s and the decimal part, i.e. n s = n° s + a, where the 
number a satisfy —1/2 < a < 1/2, which can be rewritten as, 

_ 1 10c l( 4n + l)-c 2 (2, i + 3) 1 
s 2 80ci - 4c 2 s 2 V ; 

(a) Because the singlet "pair" number isn't negative, for n° < 0, n s must be zero. Hence, 
in this case, the ground state is \n, 0, 2n, 2n\ A >. 

(b) < n° s < ~. In this case, the eigenenergy is lower when the singlet "pair" number is 
n° s than any other integer. So the ground state is \n — 2n° s1 n s , 2n — 4n°, 2n — 4n°; A >. 
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(c) For n° > ~, the singlet "pair" number takes the highest value as it can. With differ- 
ent total bosons number, there are two follow: (i)when n is even, the ground 
state is |0, n/2, 0, 0; A >; (ii)when n is odd, the ground state is |1, n/2, 2, 2; A >. 



IV. THE SUPERFLUID-MOTT INSULATOR TRANSITION 

In this section, we consider the hopping term and calculate the first and second or- 
der corrections to the ground energy. Although the eigenstates, among whose quantum 
numbers only magnetic quantum number is different, are degenerate, we can do perturba- 
tive calculations regarding the degenerate states as nondegenerate ones for the off-diagonal 
hopping matrix elements between states with the same Boson number n are zero, i.e., 
< n ,n s , F, F z ; X\hi\no,n s , F, F' z ; X >= 0, when F z ^ F' z . Therefore, the first and second 
order corrections to ground energy are expressed as follow 



e 



(i) 



-■<g\h 1 \g>=ztY t <l>>a « = -2,...,2 (21) 



£ g — 2^ (o) (o) ~ 2^ 2^ (o) (o) "~ 

nj^g £g £n n^g a Eg En 

The ground energy modified by adding the first and second order corrections becomes 

£ fl = 4 0) +4 1)+ 4 2) = 4 0) +**£4*(n,&c ,c 1 ,p 2 )## Q a = -2,. ..,2 (23) 

a 

where A a (n, ft, Co, c±, c-z) is related to the first and second order corrections of the spin component 
with magnetic quantum number a to zero-order ground energy. It depends on n, p,, Cq,Ci and 
c 2 , where p, = fi/zt,c = c /zt,Ci = Ci/zt,c 2 = c 2 /zt. Minimizing the ground energy function 
modified in second order, we find that (f> a = when A a (n, p, Co, 5i, C2) > and that (f> a ^ 
when A a (n, p, Co, ci, c 2 ) < 0. This means that A a (n, p, cq, cx, c 2 ) = signifies the boundary 
between the superfluid and the Mott insulator phases of the spin component with magnetic 
quantum number a. 

Equiping the perturbative calculations we can draw phase diagrams. The phase diagrams 
show that there is a phase transition from Mott insulator with integer number bosons to su- 
perfluid when the ratio co/t is decreased to a critical value. In the zeroth order, i.e., neglecting 
the hopping term, the ground state is Mott state in which the occupation number per site is 
pined at integer n = 1,2, corresponding to a commensurate filling of the lattice. Different 
Mott ground states maybe contain different spin components. For example, there is only spin 
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component with Zeeman level m = 2 when occupation number per site n — 1; spin components 
with m = 0, ±1,±2 for Mott state FU^o 1° >)< and 

spin components with m = 0,1,2 for 
IL(4 2) |o >)j. It is easy to realize that only one superffuid component appears when lowering 
the ratio Co/t for the initial Mott ground state containing only one spin component such as the 
case n — 1 in FigHHland n — 1, 2, 3 in Fig|SJ For the Mott ground states containing multiple 
spin components, when lowering the ratio Co/t, multiple superffuid components appear. The 
phase boundaries between superffuid and Mott insulator phase for different spin components 
are identical for some Mott ground states such as n = 2 in Fig. 1,3,4,6. and n = 3 in Fig. 1,2,4, 
and different for some Mott ground states such as n = 2 in Fig. 2 and n = 3 in Fig. 3, 6. 

We find that the position of phase boundary is related to average occupation number of spin 
component in the initial Mott ground state, i.e., the larger the average occupation number of 
spin component per site is, the easier the transition from Mott insulator to superffuid phase. 
For example, in FigQ the average occupation numbers of spin component with Zeeman level 
m = 2 are 1, 2/5 and 3/5 per site for n = 1,2,3 Mott ground states respectively, and the 
minimum critical value of Cq are 5A2zt, 2.76zt and 3.54zt for n = 1,2,3 Mott states. We can 
arrive at the same result by analyzing other figures. 



V. CONCLUSION 



In this paper, we have investigated the quantum phase transition from Mott insulator to 
superffuid phase of spin-2 cold bosons at zero temperature. First, we diagonalized the Bose- 
Hubbard Hamiltonian without hopping term and got the eigenvalue and ground states. Then 
taking the hopping term as a perturbation we calculated the first and second order corrections 
and drawn the phase diagrams, which show that there is a phase transition from Mott insulator 
with integer number bosons at each site to superffuid phase when the ratio c$/t is decreased 
to a critical value. Different Mott ground states maybe contain different spin components. 
For some Mott ground states, the different spin components appear at different moment in 
superffuid phase when lowering the ratio co/t. The position of phase boundary is related to 
average occupation number of spin component in the initial Mott ground state. 
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FIG. 1: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = O.lzt and c% = O.lzt. Here co/zt and ji/zt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superfluid phase 
and Mott insulator phase respectively. 
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FIG. 2: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = 0.02zt and C2 = O.lzi. Here co/zt and \ijzt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superfluid phase and 
Mott insulator phase respectively. For n = 2, the interior line is the phase boundary of spin component 
with Zeeman level m = 1, and the external double lines, which are too close to be distinguished in 
the diagram, are that of spin components with m = 0, 2 respectively. The inset shows an expansion 
of the part of the phase boundaries for n = 2 in the dotted frame labelled by A . 
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FIG. 3: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = 0.02zt and C2 = — 0.1 zt. Here co/zt and fx/zt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superfluid phase 
and Mott insulator phase respectively. For n = 3 case, the interior line is the phase boundary of spin 
component with Zeeman level m = 2; the middle line the phase boundary of spin component with 
Zeeman level m = —2, and the external triple lines, which are too close to be distinguished in the 
diagram, the phase boundaries of spin components with Zeeman level m = 0, ±1. The inset shows an 
expansion of the part of the phase boundaries for n = 3 in the dotted frame labelled by B. 
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FIG. 4: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = O.lzt and C2 = —O.lzt. Here co/zt and [ijzt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superfluid phase 
and Mott insulator phase respectively. 
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FIG. 5: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = —O.lzt and c% = O.lzt. Here co/zt and fi/zt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superfluid phase 
and Mott insulator phase respectively. 
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FIG. 6: Phase diagram of Bose-Hubbard Hamiltonian obtained from second-order perturbation theory 
with solid lines for c\ = — 0.02zi and C2 = —O.lzt. Here co/zt and \ijzt are dimensionless. The dashed 
lines indicate the zeroth order phase diagram. In the diagram, SF and MI denote superffuid phase 
and Mott insulator phase respectively. For n = 3 case, the interior line is the phase boundary of spin 
component with Zeeman level m = 2; the middle line the phase boundary of spin component with 
Zeeman level m = —2, and the external triple lines, which are too close to be distinguished in the 
diagram, the phase boundaries of spin components with Zeeman level m = 0, ±1. The inset shows an 
expansion of the part of the phase boundaries for n = 3 in the dotted frame labelled by C. 
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